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Abstract. Any connection between dark matter and extra dimensions can be cognizably 

bi)| evinced from the associated effective energy-momentum tensor. In order to investigate and 

test such relationship, a higher dimensional spacetime endowed with a factorizable general 

metric is regarded to derive a general expression for the stress tensor - from the Einstein- 

Hilbert action - and to elicit the effective gravitational potential. A particular construction 

^T) • for the case of six dimensions is provided, and it is forthwith revealed that the missing 

r^ . mass phenomenon may be explained, irrespective of the dark matter existence. Moreover, 

the existence of extra dimensions in the universe accrues the possibility of a straightforward 

mechanism for such explanation. A configuration which density profile coincides with the 



m 

pg . Newtonian potential for spiral galaxies is constructed, from a 4-dimensional isotropic metric 

plus extra-dimensional components. A Miyamoto-Nagai ansatz is used to solve Einstein equa- 
tions. The stable rotation curves associated to such system are computed, in full compliance 
to the observational data, without fitting techniques. The density profiles are reconstructed 
^ ', and compared to that ones obtained from the Newtonian potential. 
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1 Introduction 

General relativity (GR) provides solutions to the Einstein equations that currently offer 
important support to understand, for instance, black holes and cosmology. Mainly in the 
last two decades, such formalism was extended for extra-dimensional spacetimes. There are 
many motivations to consider such approach: from the Kaluza-Klein - and string theory as 
well - idea of unified fields, to the recent effort to deal with the hierarchy problem from 
the idea of diluting gravity in submillimetric extra dimensions [1]. Another motivation - 
indeed more fundamental for GR - is that higher dimensions may play a fundamental role 
to unravel profound interpretations concerning the nature of spacetime. On this behalf, a 
series of works deal with a (4 + n)-dimensional metric splitting, in order to find solutions for 
the Einstein equations, in both static or stationary spherical/axisymmetric geometries. The 
static, spherically symmetric solutions of Kaluza-Klein theory were calculated independently 
by several authors [2]. For instance, in [3] and [4] some stationary spherically symmetric solu- 
tions on Kaluza-Klein theory are provided, as well as in 5-dimensional GR [5]. Furthermore, 
in [6] some asymptotically flat Kaluza-Klein solutions, corresponding to regular black holes 
in four dimensions, are regarded, and prominent features concerning thermal evaporation are 
presented. Such formalism is generalized in a (4 -|- n)-dimensional scenario [7]. It is accom- 
plished in the so called abelian Kaluza-Klein theory — possessing the U(l)" as its internal 
isometry group — which lower dimensional case was investigated in [8]. The rotating dyonic 
black holes of Kaluza-Klein theory were investigated [9] withal, and the explicit connection 
between such models and extra-dimensional theories is explored in [10]. 

Therewith, one of the most outstanding non compact extra dimensions usage is the 
Randall-Sundrum (RSI) model, wherein the so called braneworld is embedded in an AdSs 
space [1]. Moreover, some large extra dimensions formalisms were investigated in a (4-|-l-|-l)- 
dimensional setting. The Einstein equations with negative bulk cosmological constant and a 
3-brane, with appropriately tuned energy momentum tensor, were solved [11]. In the case of 



two extra dimensions, which shall be considered as a particular case of our formalism, some 
aspects were approached at [12]. 

The pivotal physical motivation in our formulation is that gravity with extra dimensions 
can possibly shed new light on the so called "dark matter" puzzle [13-16]. The idea is not a 
novelty, since some theories, as universal extra dimensions [17], and some variants of Kaluza- 
Klein theories as well, present particles corresponding to extra-dimensional excitation modes 
- representing a dark matter candidate (for a review see [18]). Here the approach is different, 
in the sense that it is argued whether the existence of higher dimensional spacetimes suffices 
to explain the missing mass problem {aka "dark matter"), without any specific dark particle. 

In order to precisely ascertain what this assumption infuses in the dark matter problem, 
we delve into a generalization where our universe has D = 4 + n dimensions. Moreover, at a 
first glance nothing precludes the matter to live - or leak as well - in the n extra dimensions 
(see e.g. at the level of particle physics, the UED proposition [17]). In what follows the 
formalism is developed a la Kaluza-Klein 

Accordingly from such scenario, an Einstein-Hilbert gravitational action 

S = Jd^X d^y^-i^+n)g ((^+-^ R + Cm) , (l-l) 

where Cm denotes the matter Lagrangian, leads to the field equations 

^'+''^Gab = ^""^^^Tab, (1.2) 

where A,B = 0, 1, . . . , 4 + n, and y denotes the extra dimensions. Hereon the tensor com- 
ponents ^^'^"''Gab and ^^'^"''Tab are respectively denoted by Gab and Tab (the same for 
the curvature tensor and scalar, and for the metric components as well). No modification is 
considered for the gravitational constant G. 

Given such approach, in Section 2 the spacetime metric is split, in order to express the 
4-dimensional stress tensor. It infuses some interpretations about the extra terms, appearing 
due to the inclusion of new dimensions (Section 3). In Section 4 the equations of motion are 
developed and in Section 5 we calculate - by linearized gravity - the Newtonian limit and a 
modified Poisson equation for the theory. The Section 6 is dedicated to show how one can 
interpret dark matter as a prominent pure effect, coming from the spacetime extra dimensions. 
We provide particular cases where two extra dimensions are regarded. Stable galaxy rotation 
curves are presented together with their associated analytical expression when the system is 
axisymmetric in Section 7. Finally in Section 8 some results with real galaxies are presented 
and analyzed what, together with the results in Section 7, put our formalism in compliance 
with observational data. 

A companion paper shows a phenomenological example [16], as well as other semi- 
phenomenological ones provided by [14]. The assumptions here are straightforward, and 
quantum field theory is not explicitly taken into account. GR is the theoretical limit con- 
cerning the present work. In this article, unless explicitly mentioned, c = 1 = 87rG, and the 
signature adopted for the 4-dimensional part of the metric is diag(— , -|-, -|-, -|-). In addition, it 
is considered ab initio that the extra-dimensional part is constituted by spacelike coordinates. 



2 Splitting the metric 

Our universe is supposed to have 3 + 1 + n dimensions. Initially, the extra dimensions 
compactification is kept apart. The most general metric for such universe is given by 

{9al3 gab\ /r, 1 N 

gAB = [ ^ „ h (2-1) 

\9al3 9abJ 

where a, /3 = 0, . . . , 3 and a, 6 = 4, . . . , n. Furthermore we consider the convention to make 
the metric as a function of 3 + 1 coordinates only: qab = dAsix"')- This metric components 
QAB contain the 3 + 1 universe metric terms g^i^ and the extra dimensional terms gab, as well 
as the crossed components. Eq. (2.1) can be rewritten for convenience as 

9AB = galdSA^B + aab^A^B + gab^A^B + gafiSA^s, 

where (5'- are the Kronecker symbols. The derivatives for such metric components are given 

by 

gAB,C = gal^n^A^B^C + gab^^A^B^C + gab^^A^B^C + dap^-y^A^B^C- (2-2) 

The case qab = gap^'X^B + gab^A^B is considered here, motivated by formalisms where 
g^ = — namely Randall-Sundrum warped metrics [1] and other physically relevant cases 
[10, 13, 19, 20, 22]. In particular, a comprehensive program on observations and measurements 
on braneworld effects in astrophysics and particle physics is provided in [21]. Furthermore, 
some generalizations regarding variable tension braneworld models and cognizable physical 
effects were presented [22]. 

The inverse metric is written as 

gAB^gaP^A^B^gab^A^B^ (2.3) 

and the derivatives are straightforwardly provided by Eq. (2.2) 

gAB,C = gaP,^5'X5^B^l + QabrfSA^^B^c- (2-4) 

Assuming that the spacetime has a connection presenting no torsion, one yields the following 
Christoffel symbols {^} = \g {gBM,c + gcM,B — gBC,M)- Splitting this last expression 
by Eqs. (2.3) and (2.4) it reads 
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{B^h} = {^" } syB^l + 2 [9''"'i9bm,X^%sl + gcm,pd^€s^c) " g""' gbc,X^%sh\ ■ (2.5) 

The Ricci tensor components are well know to be written as 

Rab = du {ab] - 9b {am} + {ab} {n&i} - {am} {nb} ■ 

Taking into account that the terms of the metric depends solely on x°, the equation above 
reads 

Rab = RapS'XS^B + RabS'Xd'h ■ (2.6) 

The splitting (2.3) and the derivatives (2.4) of the metric in the Christoffel symbols (2.5) 
lead to the following expressions for Eq. (2.6) 



Rap — T^ap — -^ 



g ^pgmc,a i" g gmc,aP ~r — y gmc,ag gne,p 



(2.7) 

(2.8) 
where TZap = d^ij;'^} - d/s {a%} + {afs} {Jli}-{au} {fp} is the conventional (3+1) Ricci tensor. 



Rab - "2 



9%5a6,7 +5^'^9afe,7M " 7T {g''''gac,tMg''"' gnb,'yg''"' gam,-yg"'''gbc,,y - g""' gab,,yg"'''grnc,'y) 



3 The energy-momentum tensor 

Consider now the Einstein tensor components Gab = Rab — i^/'^)RgAB, where R = 
g^^ Rmn is the scalar curvature. For a stress tensor written as Tab = T'q,^(5^(5^ + Tab5\5'^^, 
one can follow the same arguments accomplished for Rab- Introducing the results obtained 
in the previous Section in Eq. (1.2), the split components of the stress tensor is given by: 

Tap = Rap - -^{g RMN)gal3, Tab = Rab " ^(^ RMN)gab- (3-1) 

The tensor Tap represents the energy/pressure content in the (3 + l)-dimensional landscape. 
This is the part that is clearly of paramount interest. Regardless of Tab 7^ 0, here it is 
assumed that only the Tap visible part can bring some light to observational issues. From the 
expression developed for Rap, Eq. (2.7), the influence of extra dimensions is evident. From 
(2.7), (2.8), and (3.1) it yields 

Tap = Tap + la/3, (3.2) 

where Tap represents the (3+l)-dimensional influence, and 1-ap denotes the correction elicited 
uniquely from extra dimensions, where explicitly 

Tap = IZap - \{g'"'n^u)gap, (3.3) 



la/3 = -^ \^ap - ^ig'"''n^u)gap + (5"'"^mn)5a/3 



(3.4) 



Note that the stress Tap has the conventional form for the Einstein tensor Gap- The tensor 
yiap represents a "curvature" term evinced from the presence of extra dimensions, and is 
given by 

^ap = g^%gmc,a + g^^ gmc,ap + -^g^^ gmc,ag^''^ gne,p- (3.5) 

In addition, from the action (1.1) the stress tensor can be derived from the conventional 
definition Tab '-= —2 ^ ^^ + gAB^M- In this new approach, it is clear that the (3 + 1)- 
dimensional part Tap is proportional to the Lagrangian multiplied by gap and some quantity 
from the variation bjbg^^. 

4 An expression for the equations of motion 

A general Lagrangian for the gravitating test particles in a spacetime with extra dimen- 
sions, with metric elements gAB = gap^'X^B + gah^\^\, can be written as 



L = [gABi^i^'f^ = {gapx'^xf' + gabx'^x'Y'^ (4.1) 

re X = dx /ds. Tl 
^7 = 0. As Oa = dab A + da5\ and gab = gabix"") it follows that 



where x = dx /ds. The motion equations come from the Euler-Lagrange expression -^ {'^^ 



dL _ dL ^ dL 
dL 



dx"i 2 



-L ^{gap,'yX°'x^ + gab^x'^X^), 



and ^ = 0. It immediately yields 



9xC 2 

Likewise, the term ^ {'§±c) can be developed: 

dL 1 



dL 1 






Now 



d r dL 



ds \ dx"i 



L 



-1 



dx° 



•-^ ] X- x^ + g^^x' 



Since x"" are cyclic variables, it reads the integration constant: 

g.mx'^ = N,, (4.2) 

a constant vector. Hence ^ (^p) = 0- Inserting the terms together, multiplying by Lg'^'^ and 
using (4.2) the equations of motion are derived: 

^^ + {a^} i"^^ = \9ab,,9''^N,g'''Ndg'''. (4.3) 

Clearly it is possible to realize that the extra dimensions induce an external force in the 
system, that depends on gab and Nc- 

5 Equation for the visible field 

The main aim now is to compute the gravitational potential in the Newtonian limit, 
since galaxies and clusters can be described physically as Newtonian objects — corresponding 
to the approximation in which gravity is weak. The weak limit is assumed uniquely in the 
4-dimensional spacetime: the deviation ^a/s of the 4-dimensional metric ga/s = rjajj + 70,^3 
is small {rjap denotes the Minkowski metric). For the linearized gravity, the stress (3.3) is 

essentially written from linearized curvature Tl)y.a = ^^ 1 afl f ~ ^^ 1 uS f ^^ 

TalS = -^9^5^7«/3 + 9^9(^7a)M - l^afsd'^d'^^.u (5-1) 

where 7q,^ = 70-/3 — ^rjapj is the traceless part of ^afs and 7 = j^. Linearized gravity has a 
gauge freedom given by 7^/3 1— )■ 7o/3 + i^g??a/3, where £^ denotes the Lie derivative with respect 
to the generators ^^ of a differential diffeomorphism. To the first order, such transformation 
represents the same physical transformation as 70,^3 . This gauge freedom is used to simplify 
the linearized Einstein equation. Solving the equation d^dp^a. = —d^lap for ^q, a gauge 
transformation [23] that leads to d^j^^a = — similar to the Lorentz gauge condition — can 
be elicited to obtain the simplified Einstein equation 

Tai3 = - j9^9M7a/3- (5.2) 
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Figure 1. (a) The rotation curves of an effective potential where the Poisson equation is solved from 
a Miyanioto-Nagai ansatz, with parameters a = 1 and b = 0.1, describing a disk galaxy, and where 
the Laplace equation for the extra field is solved by a Chazy-Curzon disk, with stable cut parameter 
c = 1.5. They correspond to gravitating particles which orbit in a background. The full line is the 
curve when there is no extra dimensions, i. e., TV^j^ = Nz^ = 0. The dotted lines represent a set of 
extra dimensional solutions for fixed Nz-^ = 0.1 and Nz^ varying from 0.8 to 0.95, as in [14]. The 
parameters a, 6, c arc chosen in order to the orbital motion be stable, (b) The same set of solutions 
presented in (a), but now we represent the squared epicyclic frequency k^. When this last is positive, 
we guarantee the stability of the system. 



For the extra part it reads 

'ZaB 



-,{g'^''d^^d^g^n)9alS-9"'"9mn,a^ 



(5.3) 



When gravity is weak, the linear approximation to GR should be valid. There exists a global 
inertial coordinate system of r]^^ such that 



fa/3 



Tal3 + Ta/3 ~ ptatp. 



\d'd,^^p + \ 



-(ff™"a^a^5™,„)^a/3-5"'"5. 



'nin,a[S 



Ptats, 



(5.4) 
(5.5) 



where t^ is the time direction associated with this coordinate system. This equation can be 
interpreted as the modified Poisson equation considering a universe with more than 3 + 1 
dimensions. 

6 Extra dimensions and the interpretation of dark matter 

Define j^a = —4:(f), where (j) = 4'i^) is a 3-space scalar field. Furthermore consider a line 
element ds^ = Yli=i e'^'<izf, where ds^ is the world line for the extra sector, Zj denotes the 
extra coordinates and ipi = ipi{x) are potentials associated to extra dimensions. 

If one asserts, as a first approximation, that Rab = (that is not similar to impose 
vacuum in the extra space, since now Tab becomes Tab = ~\{9^^ RMN)gab)^ Eq. (2.8) implies 
at first order the sigma model g^'^ {cr ^^a~^) ^y = for the extra part, where a denotes the 
diagonal matrix representing the metric associated to the system, which implies 



d^'d^gab = 0, 



.1 



yielding the following equation 



Idf^d^^^p - ^g"'''gmn,ap = pt^tfs, (6.2) 



or in other words 



V^0 = p. (6.3) 



It means that our visible matter density profile is provided uniquely by the 4-dimensional 
field. If extra dimensions have some impact on the theory, certainly it is not as extra matter, 
but it appears at least as the term \g^"' gmn,a0 in (6.2). In matter models applicable to a 
galactic context pressure arises from the velocity dispersion in the motion of particles or from 
exchange of momentum among particles through collisions — ionized interstellar baryonic 
gases. An extra pressure for the same energy-density could be possible for the star gases, 
which have no equation of state, but would introduce modifications on the equation of state 
of the interstellar gases. Instead, the extra term in equation (6.2) has the form of a tidal force 
produced by the purely geometric effect of the extra dimensions. Likewise, the equation of 
motion does not necessarily dismiss the influence of the extra field, as it can be illustrated 
by the following example. A general metric can be written for the simplified case of six 
dimensions as ds^ = — (1 — 2(/))di^ + dx.dx + e~^dzf + e^dz^ [14], where dx • dx is the 3- 
dimensional line element. Now, in the Newtonian limit the motion is conceived to be much 
slower than the speed of light, and i;° can be assumed to be (1, 0, 0, 0) in (4.3), as well as the 
affine parameter s may be approximated by the coordinate time t. Thus it follows that 



d^r^ 1 

"•" too/ 



dt2 ' LUOJ 2 



d{e-^) 



,2V Ar2 



K + 






-^^Nl \ . (6.4) 



As time derivatives of (f) and V' are neglected, Eq. (6.4) leads to 

a = -V$, (6.5) 

^ = c^+\{e-^Nl+e^Nl), (6.6) 

that is the equation of motion of test bodies gravitating in an orbit in a given background for 
our system with an effective gravitational potential $. To find a form for those potentials, 
from (6.1) and (6.3) it yields 

V^V -Vip-Vip = Q (6.7) 

V^ = p. (6.8) 

Non-linear terms do not appear, since the a matrix is diagonal. In particular, Eq. (6.7) can 
be rewritten as 

V'x = 0, (6.9) 

where the identification x = ^ is accomplished, and hence Eq. (6.6) has the form 

<^ = <t^ + \{xNl + X-^Nl). (6.10) 

Together with Eqs. (6.8) and (6.9), it represents a complete system of equations for a gravi- 
tational effective potential in 4-dimensions, assuming a 6-dimensional universe. Note as well 
that (6.9) is the Laplace equation in flat 3-dimensional space, and therefore x can be taken 
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Figure 2. Values for the general 6-dimensional model (Newtonian limit): N^^ = 0.2, Nz^ ~ 0.8 and 
0.9, a = 1, b = 0.01. The violet lines are two stable examples of the present model. The full blue line 
is the relativistic model presented in [16]. 
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Figure 3. Values for the general 6-dimensional model: N^^ = 0.18, N^^ = 0.85 and 0.9, a = 1, 
b = 0.01. The violet lines are two stable examples of the present model. The full blue line is the 
relativistic model presented in [16]. 



as a solution of Laplace equation for an appropriate Newtonian source of any symmetry. 
The form of (j) can be indeed expressed as the point source potential ~ 1/r for a Newtonian 
system. 



In a system which acceleration a, Eq. (6.5), has radial direction, e.g. an idealized system 
similar to a galaxy, the rotation curves for a potential 4> ~ — 1/r are provided by 



Vc 



- + -rtl;^,{Nle- 



Nle^) 



.1/2 



(6.11) 



7 Solutions 

Now we look for a way to implement an exact Newtonian solution for the Poisson 
equation (6.8), as well as for the Laplace equation (6.9). As a system similar to a galaxy is 
regarded, the Poisson equation can be implemented e. g. by a Miyamoto-Nagai potential, and 
for the case of Laplace equation it can be used a Chazy-Curzon solution for an axisymmetric 
configuration. Namely the Miyamoto-Nagai ansatz is written as [24] 



^{R,z) 



M 



R^ + {a + Vz^ + fo2)2 



(7.1) 



where o, b are positive constants. The Chazy-Curzon solution for a particle of mass m in the 

position z = zq is given by [25, 26] 

2?ra ,„ , 

X = ^, (7.2) 



where R = ^r^ + (z — zq)'^. In this last, it is possible to fix zq = 0, and to introduce the cut 
method to generate a disk solution [27], such that z i— )• |2;| + c, where c is the cut parameter. 
Fixing values for a, b and c in such a way that b/a ~ 0.1 - providing a similar light distribution 
of a disk galaxy [27] ~ and c > 1 (by stability issues and to prevent relativistic disks [14]) it 
is possible to find the set of curves given in Fig. la. 

The stability of this last configuration is guaranteed by the positive epicyclic Newtonian 
parameter k^ = d'^^/dr'^ + SVj/r^, as shown in Fig. (lb). 

A more useful manner to analyze the effective potential is writing it as 



$ = (f> + Ccosh{ilj + 6), 



(7.3) 



where C and 5 are constants to be determined. This is the correspondent Newtonian potential 
in the present formalism. Clearly, from Eq. (6.6) those constants are related to A^^^ and Nz2 
as 



C 



Ar2 _7v2 

Z2 Zl 



tanh 5 = —-5- . 

A£ 



Given the present results, one can ask about the nature of the potential (7.3). The term 
comes exclusively from extra dimensions and it may be related to a "dark matter" potential 
as 

(/'extraD = <pDM = CcOsh(^ + 5), (7.4) 

showing that our approach is equivalent to a dark matter effect, although the concept above 
is completely apart from the concept of extra matter. 



8 Probing results with real galaxies 

Taking into account only the stable curves, in Figs. 2 and 3 we compare with some 
optically observed rotation curves of spiral galaxies. 

Here we are not providing a composition of halo dark matter velocities plus the disk 
gas and the velocities of stars. What is happening is that the clean stable calculated curves 
are simply fitting the region of interest. The not surprising ad hoc adjustment of Nz-^ and 
Nz2 actually could assert nothing about the astrophysical role of the extra dimensions in the 
model. However, the calculation of stable configurations brings realizable values for A'^^^ and 
Nz2, which makes it possible to visualize a minimum representation of a real disk galaxy. 
We furthermore compare to some phenomenological models used in astrophysics as Navarro, 
Frenk &: White [28]. The observed data is taken from [29-31], but for an alternative data 
source see [32] , where our model is also successful. 

Note that our results reproduce with great fidelity the shape of non-planar curves, as 
that appearing in Fig. 3(b) (data taken from [33]). 

9 Concluding remarks 

At the present work we showed how to calculate a "dark" potential from extra-dimensional 
imprints. If extra dimensions have some impact on the theory, certainly it is not as extra 
matter, but it appears at least as tidal force produced by the purely geometric effect of the 
extra dimensions, namely the 2ff™"5mn.,a/3 term in (6.2). This is confirmed by the extra force 
appearing in the equations of motion. 

The Newtonian limit is calculated and we accomplished the visible potential that has 
two main components: a term from 4-dimensional metric and another term coming exclusively 
from extra-dimensional effects. The 4-dimensional term obeys the Poisson equation and the 
extra term is solved from a Laplace equation. Two situations are provided: a first approxi- 
mated model and a second that arises from exact solutions of Laplace and Poisson equation, 
respectively the Chazy-Curzon disk solution and Miyamoto-Nagai ansatz. Those examples 
are calculated for a case where two extra dimensions are taken into account. Moreover, any 
analogous weak field potential decomposition, as the effective ones for elliptic galaxies or 
globular clusters, in which (p would be a different potential with the appropriate symmetry, 
holds for this formalism. 

The extra field arises naturally from metric extra terms, and the modification of gravity 
shown from the calculation of rotation curves, can present an explanation about the missing 
mass effect in astronomy. The effective potential calculated here is a source to explain the 
phenomenological results found in [16] and [14]. The extra part of this effective potential 
has the form C cosh[ip{x°') + 6]. The gradient of such potential represents the extra force 
associated to the presence of extra dimensions. It modifies the rotation curve profiles and 
rise anomalies in potential associated to clusters. The Laplacian corresponds to a "dark"' 
density profile, that is negligible compared to the visible density V^^ = p (see Section 6). 
Furthermore, the extra-dimensional proposal is consistent with a a more general length scale 
and galactic morphology. 

The present calculation is an alternative approach to understand the dark matter prob- 
lem. In fact, to solve the dark matter problem means to solve a long list of cosmological 
anomalies. In this aspect, as a complement, considerations about cluster lensing are also 
addressed in [13]. 
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